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Abstract 

We discuss renormalizability of a recently established, massive gravity theory with 
particular higher derivative terms in three space-time dimensions. It is shown that this 
massive gravity is certainly renormalizable as well as unitary, so it gives us a physically 
interesting toy model of perturbative quantum gravity in three dimensions. 
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1 Introduction 



By now it is famous that it is very difficult to construct a quantum theory of Einstein's general 
relativity in four dimensions because of nonrenormalizability. The nonrenormalizability is 
traceable to the fact that its perturbative expansion parameter is the Einstein's gravitational 
constant k, with dimension of length squared. Thus, the divergences encountered in the 
perturbation theory of Einstein's gravity are totally out of control and we have no physically 
meaningful prediction on the basis of perturbative quantum gravity 

There is, however, an obvious renormalizable generalization of Einstein's general relativity 
in four dimensions, which is the 'higher derivative gravity' theory with higher derivative 
curvature terms [1], but this theory is known to be non- unitary owing to the presence of the 
non-unitary massive ghost, unfortunately. 

With such a situation, it has been anticipated for a long time that this problem might be 
cured by going to lower space-time dimensions, but it is remarkable to recall that we have not 
yet had no grasp of perturbative quantum gravity, which is dynamical, unitary and possibly 
power-counting renormalizable, even in the lower dimensions. 

To date, there is only one candidate as known exceptions to this picture, which is, what we 
call, topologically massive gravity in three dimensions [2, 3]. This geometrical gravitational 
theory is power-counting renormalizable despite being dynamical, unitary and including a 
dimensional coupling constant like Einstein's gravity in four dimensions. However, even this 
hopeful theory has a flaw in the proof of renormalizability in the sense that there is no gauge- 
invariant regularization in such a way to preserve the desirable power- counting behavior. 
Related to this fact, note that dimensional regularization cannot work well because of the 
explicit Levi-Civita tensor density e^'^^ , gauge-invariant, higher derivative regulators spoil the 
arguments of formal renormalizability, non-covariant cutoffs cannot be easily analyzed, and 
non-local regularization method involves some assumption to be proved [4, 5] . 

Recently, in three space-time dimensions there has been an interesting progress for ob- 
taining a sensible interacting massive gravity theory [6] ^. This model has been shown to be 
equivalent to the Pauli-Fierz massive gravity [16] at the linearized approximation level and 
thus massive modes of helicities ±2 are physical propagating modes. A key idea in this model 
is that one adds higher derivative curvature terms to the Einstein-Hilbert action with the 
wrong sign in such a way that the trace part of the stress-energy tensor associated with those 
higher derivative terms is proportional to the original higher derivative Lagrangian. With this 
idea, it turns out that the scalar mode coming from higher derivative Lagrangian is precisely 
cancelled out [17]. More recently, this new massive gravity model in three dimensions has 
been studied from various viewpoints such as the unitarity and the impossibility of general- 
ization to higher dimensions [18], relation to the Pauli-Fierz mass term [17], the AdS black 
hole solutions [19], the properties of linearized gravitational excitations in asymptotically AdS 
space-time [20, 21, 22] and AdS waves [23]. 

The aim of this article is to show that the new massive gravity theory in three dimensions 

^See the references [7, 8, 9, 10, 11, 12, 13, 14, 15] for alternative massive gravity models. 
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[6] is indeed renormalizable within the framework of perturbation theory. It is worhtwhile to 
notice here that we can make use of dimensional regularization in the new massive gravity 
theory since there is no the Levi-Civita tensor density e'^^'' in the action, which should be 
constrasted to the case of the topologically massive gravity. 

The reason why the new massive gravity is renormalizable is very simple. With a suitable 
choice of the regulator and gauge condition for diffeomorphisms, the graviton propagator falls 
off like ^, implying that the gravitational field has the ultraviolet dimension of (mass)" 2. 
Then, the power counting argument together with the Slavnov- Taylor identity and renormal- 
ization equation for the effective action shows that the divergent part of the effective action 
consists of terms of dimension 3 at most and is invariant under diffeomorphisms. Thus, 
only the admissible counter-terms are the cosmological term and the Einstein-Hilbert action, 
and thereby the theory is perturbatively renormalizable. Note that the higher derivative 
curvature-squared terms receive no corrections so that the property of the unitarity in the 
new massive gravity is not violated even after the renormalization procedure. 

Since it has been already shown that this massive gravity is unitary, we have for the first 
time a physically plausible quantum theory of gravity within the framework of perturbation 
theory such that it is dynamical, renormalizable and unitary although this theory can be 
formulated only in three dimensions. 

In the next section, we perform the BRST quantization of the new massive gravity theory 
in three dimensions. In the third section, we derive the propagator for the gravitational 
field on the basis of the gauge fixed, BRST-invariant action obtained in the section 2. In 
the fourth section, we consider the power-counting and determine the divergent part of the 
effective action by using the Slavnov- Taylor identity and the renormalization equation. The 
final section is devoted to conclusion and discussions. 



2 BRST quantization of the new massive gravity theory 

We start with BRST quantization of the new massive gravity theory in three dimensions 
[6] which is an interactive and unitary theory with higher derivative terms for the massive 
graviton. The classical action which we consider is of form ^: 

where = 16nG {G is the 3-dimensional Newton's constant) and M is a constant of mass 
dimension. Let us note that k has dimension of {mass)~2 ^ so the theory defined by the action 
(1) might at first sight appear to be unrenormalizable. As usual let us write 

fl'/xi^ = Vfiiy + (2) 
"^The space-time indices /i, i^, • • • run over 0, 1, 2. We take the metric signature (— , +, +) and follow the 
notation and conventions of the textbook of MTW [24] . 
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where a flat Minkowski background rj^i, has the diagonal element (—1,1,1). With this def- 
inition, the gravitational field h^i, has canonical dimension of {mass) 2. Consequently, the 
Einstein-Hilbert action and the graviton-matter interaction terms have canonical dimensions 
greater than three, which is the origin of unrenormalizability of Einstein's general relativity 
without the higher derivative terms. Later, we will see that the gravitational field /i^i^ has 
ultraviolet dimension of (mass)" 2, whose fact leads to renormalizability of the new massive 
gravity theory under consideration. 

Now the BRST transformations for the metric tensor, ghost, antighost and Nakanishi- 
Lautrup auxiliary field are respectively given by 

SbQuu = -«^(V^Ct, + V,.c^), 

Ssb^ = 0, (3) 

where the covariant derivative is defined as usual by V^Ci, = d^Ci, — T^^^ca with the affine 
connection T^^. It is straightforward to prove that the BRST transformations (3) are off-shell 
nilpotcnt. 

Next, we wish to fix the gauge invariance by some suitable gauge fixing condition. The 
usual gauge fixing condition would be the de Bonder's gauge (or its higher derivative gener- 
alizations) because it is known that the renormalization procedure becomes extremely simple 
in this gauge as first pointed out by Stelle [1]. However, it turns out that only the de Donder's 
gauge does not work well in the present theory owing to the scalar mode whose propagator 
does behave as ^ unlike A for large momenta. 

It is clear that this scalar mode is nothing but the conformal mode stemming from the 
Einstein-Hilbert action since the conformal mode from the higher derivative curvature terms 
is exactly cancelled out by a specific combination of them. The similar situation has been 
already happened in the case of the topologically massive gravity theory where there exists 
the conformal mode only from the Einstein-Hilbert action since the gravitational Chern- 
Simons term, which is of higher- derivative nature, is conformally invariant so that there is 
no conformal mode in this term [4]. Indeed, the existence of the scalar mode has also caused 
some trouble in the proof of renormalizability in the topologically massive gravity theory. 
However, at the linearized level the new massive gravity is shown to be equivalent to the 
Pauli-Fierz massive gravity with only the spin ±2 massive graviton modes without the scalar 
mode by path integral ^, so the troublesome scalar mode might be somewhat an artificial 
degree of freedom. 

In order to resolve this problem associated with the propagator for the scalar mode, we 
shall obey the following line of the argument. First, the key observation is to recall that 

^But there has not yet been the proof of this statement by canonical formahsm. 
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provided that we add the PauU-Fierz mass term 



2 

CpF = -^{h,.h^''-h'), (4) 

to the action (1), we have no propagator for the scalar mode while the propagator for the 
massive gravitational modes of helicities ±2 falls off like ^ for large momenta [17] as desired. 
However, if we simply added the Pauli-Fierz mass term to the action (1) without care, the 
BRST invariance would be broken by the additional term and as a result the unitarity would 
be violated whose situation we have to avoid from physical fundamental principle. Moreover, 
it is obvious that addition of the the Pauli-Fierz mass term would modify the physical content 
of the original theory. 

Accordingly, we shall add the Pauli-Fierz mass term to the action (1) as some (infrared) 
regulator (which means that we take the vanishing mass limit after renormalization) in a 
BRST-invariant manner. Such a BRST procedure has been already studied from a different 
interest [25] , but the method turns out to be applicable to the present problem as well. Using 
the BRST procedure in Ref. [25], it is easy to obtain the following gauge fixed, BRST-invariant 
action 

S = I d^xC 

- (d''h^,-di,h){mA^-d''ip)-ic(Dc-mdi,c^)], (5) 



/ 



where we have defined g^"" = y^g^'', D^'' = + g^^^d^d^ - g^"'^p - {^p~g^"') and F^^ 

dp,Ai, — d^Afj,. Furthermore, in order to get the BRST symmetric theory, the Steuckelberg 
fields A^^ and 99, and the new ghost c, the corresponding antighost c and Nakanishi-Lautrup 
field b are introduced. In addition to Eq. (3), the BRST transformations for them are given 

by 

SbA^" = K'^i-mc^' + d^'c), 

Sbc — ib, 

SBb = Sbc^O, (6) 

where note in particular that we can set 6%A^^ = since 5bA^ = K^md'd^c^ belongs to the 
higher-order in k. 

At this stage, we should comment on the meaning of the above action (5) in detail. As the 
gauge fixing conditions of diffeomorphisms and Steuckelberg-like gauge symmetry, we have 
selected the de Bonder's gauge d^g^'^ = and the Landau-like gauge d^A'^ = 0. Then, the 
terms containing the fields Ai^,ip,b,c and c are only quadratic in fields and therefore there 
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are no interaction terms among them. In obtaining the propagators, off-diagonal pieces such 
as < b/^h'^'' >, < bA^^ > and < h^^AP > are therefore irrelevant as 6^, 6 and A^^ never appear 
in vertices. 

As seen shortly in the next section, it will turn out that the propagator for the scalar 
(conformal) mode does not exist in the action (5). Then, one asks ourselves where one 
dynamical degree of freedom associated with the scalar mode has gone away since the number 
of dynamical degrees of freedom never changes in perturbation theory. The answer lies in the 
fact that there appears the propagator for the Steuckelberg field A^, which has one dynamical 
degree of freedom owing to gauge invariance. In other words, because of the BRST-invariant 
regulator the massless pole ^ of the scalar mode is transferred to that of the Steuckelberg field 
A^. This fact could be understood more clearly by taking the more general gauge condition 
d^A^ -fh = 0. 

Finally, in the limit of m — 0, the action (5) reduces to the gauge fixed and BRST- 
invariant theory of the original new massive gravity up to irrelevant non-interacting terms. 
This limit must be taken after the whole renormalization procedure is completed. 



3 Graviton propagator 

On the basis of the gauge fixed, BRST-invariant action (5), we wish to derive the propagator 
for the gravitational field /i^j,. To this end, it is useful to take account of the spin projection op- 
erators in 3 space-time dimensions [1]. A set of the spin operators P^'^\ P^^\ P^^'^\ p(fi,w) ^ p{o,sw) 
and p(°'"'*) form a complete set in the space of second rank symmetric tensors and are defined 
as 



P(2) 
fiv,pa 


— l^i^ p.p(^vcT + d p.a()vp) 


^B B 

c^'p.v^ pai 




p(l) 


~ lj{^IJ.p^v(T ^Ha^vp 


+ dvp'-^ixa + ^ua^^iJ.p) , 




p(0,^) 
p,u,pa 

p{0,w) 
IJ,v,pa 

p(0,sw) 
IJ,u,pcr 

piO,ws} 
IJ,v,pa 


- 9 

— c^p.v"pai 

^ B n 
^j^pv'^pai 

- ^ u B 

— ^j^iiu^pc- 




(7) 



Here the transverse operator B^i, and the longitudinal operator ujjj,i, are defined as 

1 

Vp.f '^d^di/ Tjfj^i/ u}p,i/, 
uJiuu = ^d^d^. (8) 
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It is straightforward to show that the spin projection operators satisfy the orthogonahty 
relations 

p{i,a) pU,b) _ ijij Sjab p{i,a) 



p(i,ab) p(j,cd) _ njd)cp{i,a) 
Hfjpcr fXTjXT ixv,\ti 

p{i,a) pU'bc) _ nj cab p{i,ac) 
p,V,pa p(T,\T ^ ^ pv,\Ti 

p(i,ab)p{j,c) _ cij cbc p(i,ac) /qx 

with i,j = 0, 1, 2 and a, b,c,d = s,w and the tensorial relation 

[p(2) + + p(o,.) + p(o,^)l^^^^ = ^(^^^^^^ + rj^^rj.^). (10) 

It is then straightforward to extract the quadratic fluctuations in h^i, from each term of 
the action (5) and express it in terms of the spin projection operators: 

Cr2 = 

2 
2 

= -!^/iM-[p(2) + p(l) _ piO,s) _ ^^p{0,sw) ^ p(0>"'-))]^^^^^/,P<^. (11) 

Note that a salient feature of the specific combination of the higher derivative terms is the 
disappearance of the spin projection operator corresponding to the spin scalar mode: 

4 

Using the relations (11) and (12), the quadratic part in h^y of the action (5) is expressed 
in term of the spin projection operators like 



J d'x]^h^^Vp.,pah^% (13) 



where V^j,i,^pcr is defined as 



V,.,p. = [l(2_tf_n_^2^p(2)_!^p(i) + l(a + ^2)p(o.) 



2 

m 



+ ^(P(°'^") + P^°'"^^)]M.,pa- (14) 
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Then, the propagator for h^i, is defined by 

< 0|T( V(^) V(y))|0 >= iV-^lpa^^'Xx - y), (15) 
where using the relations (9) and (10), the inverse of the operator V is easily calculated as 

-p-i ^ r ^ P(2) _ A P(i) _ ° + pio,y^) 

This expression of the graviton propagator brings us a few important information. First, 
let us recall that the spin 1 component projected by P^^^ and the spin ones by P^^''^\ p{o,sw)^ 
and p(°'^^) can be gauged away, whereas the spin 2 component projected by P*-^-* and the spin 
one projected by are physically relevant. One notable feature of the propagator (16) is 
that there is no spin component projected by p(°'*), which is the reason why we have taken 
the Pauli-Fierz mass term as the regulator of the conformal mode. Second, the propagator of 
the spin 2 massive graviton modes of hclicities ±2 falls off like ^ for large momenta. 

Finally, let us notice that as the result of addition of the Pauli-Fierz mass term, there 
appear two massive poles in the sector of spin 2 graviton, which is of form 

= -^(□-a,+)(n-u-_), (17) 



lit /l I 1( ™ )'^ 

where uj± = — ^^^^ — ^^M^, which are real numbers such that a;+ > 0,a;_ < [17]. As 
examined in Ref. [17], the presence of the negative pole a;_ might seem to induce violation of 
both unitarity and causality. But such a pathology now does not occur since in the m 
limit, the positive one of the two massive poles, uj+ remains a pole of massive graviton while 
the other negative pole a;_ changes to a harmless pole of massless graviton. In this respect, 
it is worth mentioning that the original massive gravity theory has been already shown to be 
unitary and causal in the de Bonder's gauge [18]. 



4 Structure of the divergences 

We now turn our attention to the analysis of structure of the divergences of the effective 
action. To this aim, let us first add sources Kf^^, (anti-commuting, ghost number = — 1, 
dimension = |), (commuting, ghost number = —2, dimension = 1), (anti-commuting, 
ghost number = — 1, dimension = 1) and (anti-commuting, ghost number = — 1, dimension 
= 1) for the BRST transformations of g^'^, c^, and ip to the action (5), respectively ^: 

S = f d^xt 



^Here for simplicity we have regarded g^" as a basic gravitational field. 
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J Ki tv tv 

77?^ 1 1 

- -^{K^h*^^ - - -Fl, - ^hd,A>^ - (a^V - dMrnA'^ - d^v) 

1 771 

- ic{nc-md^c^) - L.c'^d^c'' + -MJ-mc'' + d^'c) + —Nc]. (18) 

K K 

Next, based on this action, let us consider the superficial degree of divergence for IPI 
(one particle irreducible) Feynman diagrams. Then, it is convenient to introduce the fol- 
lowing notation: = number of graviton vertices with two derivatives, tir2 = number of 
graviton vertices with four derivatives, = number of ghost vertices, Uk = number of K- 
graviton-ghost vertices, ul = number of L-ghost-ghost vertices, Iq = number of internal ghost 
propagators and Ie — number of internal graviton propagators. Note that the fields A^, (f, b, c 
and c are free so we can exclude such the fields from the counting of the superficial degree of 
divergence. Using this fact and the above notation, the superficial degree of divergence for 
an arbitrary Feynman diagram 7 can be easily calculated to be 

a;(7) = ^n,d, + (3-4)7B + (3-2)7G-3(^n,-l) 

= 3-nR + nR2 - no - 2nK - Shl - Ie + Ig, (19) 

where di denotes the number of derivatives in the interaction terms. Here we have made use 
of the fact that graviton propagator behaves like for large momenta as mentioned in the 
previous section. 

Furthermore, using the relation 

2nG + 2nL + uk ^ 2Ig + E^ + E-^, (20) 

with the notation that Ec — number of external ghosts and Ec — number of external 
antighosts c^, Eq. (19) is cast to the form 

3 11 

u;(7) = 3 - Tiij - (Ie - nR2) - -uk - ul - -E^ - -E^. (21) 

However, as noticed before by Stelle [1] , the graviton propagator vanishes when multiplied 
by the momenta p'*, so the terms with the de Bonder's gauge term d^g^^'' do not connect with 
the graviton propagator, and thus we can neglect such the terms. Consequently, since the 
dimension of the ghost and antighost is increased by 2, instead of (21), the correct superficial 
degree of divergence is given by 

3 5 5 

uj{^) ^S-ur- {lE-nR2) - -riK -riL- -E^ - -E^. (22) 

The superficial degree of divergence (22) gives rise to some useful information. First, since 
there is a relation Ie > uri for IPI diagrams 7, we have 

a;(7) < 3, (23) 
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which indicates that the bound on the superficial degree of divergence for IPI diagrams 7 is 
cubic to all orders, and the theory is power-counting renormalizable. In this context, it is 
worthwhile to point out that we can make use of the dimensional regularization in evaluating 
various amplitudes if necessary. 

Second, we can understand a convergent property of some IPI Feynman digrams via Eq. 
(22). In order to describe this fact in the physical terminology, it is useful to define the 
effective action F of this theory in a conventional way and expand it in a power series of the 
loop number like 

F = F(°) + F(i) + F(2) + . . . , (24) 

where F*^") describes the n-loop part of the effective action. Moreover, we separate F^") into 
the divergent and finite parts like 

r(") = rS + rg. (25) 

Now notice that all the IPI diagrams involving external ghosts or K- vertices or L-ones are 
finite since, for instance, if uk = 1, we have Ec — 1, and thus Eq. (22) leads to a; (7) < — 1. 
Hence, we can describe these convergent properties in terms of the divergent part of F^") as 

rp(n) rTn(n) c-p{n) 

div _ div _ div _ g /2g-\ 

Next we move on to the Slavnov- Taylor identity. The BRST invariance of the action (18) 
enables us to derive the Slavnov- Taylor identity for E as follows: 

= (^bE 

" 6g'"' Sc"' Scr 

9/ „ r^,, N <^E 9 5E ., 5E , 
+ K^(-mc^ + d^'c)—— + K^mc— + ib—. (27) 
^ ^5Ai^ S(fi Sc ^ ^ 



With the help of relations 



SK^ ' SL._ 

1 , „ ^„ , 5E m 
— = -(-mc^ + a^c). - = -c. (28) 



Eq. (27) can be rewritten as 

= 5rE 



o (5E (5E , 5E 5E , SE 
= i^T7:^Tz ^'^Ti— T hz6ri^ 

^ 6t St 6t St 

SMf.SAi' 5N Sip Sc ^ ' 
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If we define E by 



E = E - y d^xi^r^d^b, - ^bd^A>^], (30) 



we can obtain the Slavnov- Taylor identity for E: 



5E SE ^5E _5^5E_ SE5E_ 



where we have made use of the equations of motion to the antighosts: 

^ (5E (5E 



(5E 




(5E 




+ 












SE 




+ 


IE 



(32) 

Because at the zero-loop order we have a relation r(°) = E, it is expected that the ST- 
identity for the effective action takes the same form as (31) where E is now replaced with F. 
In fact, this statement can be easily verified. Thus we arrive at the the ST- identity for the 
effective action 

+ 7i~77;; + tttttt;: + 77^T7 ~ 0. (33) 



Let us recall that the ST-identity for the effective action naturally leads to the renormal- 
ization equation to the the divergent part F*^") of the effective action 

= 0, (34) 

where the operator Q is defined as 

5F(o) S 5F(o) S 5F(o) S 5F(o) S 



5I^_5_ 5I^_5_ 51^^ (5F(°) 5 
^ J^SK^^~5^SlZ^lAi^SM~^^~S^SN' ^^^^ 



Using the relations (26) and the obvious identities 



X-pin) rp(n) rp(n) rp(n) 

div _ div _ div _ div _ g ^gg-j 



5Ai^ 5(p SM^ 5N 
the renormalization equation becomes 



= ^?fS = TT^-r^ = D^^^cP-^. (37) 
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This renormalization equation clearly shows that F^"^ is a BRST-invariant functional of g^'^. 
Since the divergent part is local and of dimension 3 at most, the only possible form of T^^l 
reads 

= / d'x[a^n)V^ + hn)V^R]. (38) 

These divergences can be absorbed by renormalizing the Newton's constant and adding the 
cosmological constant's counter-term. Of course, after renormalization, we take the limit 
m ^ and recover the new massive gravity theory up to irrelevant free terms. 

In this way, we have completed the proof of rcnormalizability of the new massive gravity 
theory in three dimensions. One important remark is that the higher derivative curvature- 
squared terms receive no quantum corrections since they have dimension of {massY, thereby 
making it possible for the new massive gravity to make sense of even after renormalization. 



5 Discussions 

In this article, we have presented a proof of rcnormalizability of a recently proposed new 
massive gravity theory in three dimensions [6] . Our proof relies on the existence of the BRST 
symmetry and the BRST-invariant infrared regulator. 

Since it has been already shown that the new massive gravity theory is an interactive and 
unitary theory for the massive gravitons of hclicities ±2, our proof insists that this theory is 
in addition renormalizable in the perturbation theory. Thus, we have obtained a nontrivial 
toy model of perturbativc quantum gravity. 

A peculiar feature of our proof is that we have adopted the Pauli-Fierz mass term as the 
(infrared) regulator. Then, a natural question arises whether or not this method could be 
applied to the proof of renormalizability of the topologically massive gravity. We wish to 
clarify this problem in a future publication. 
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